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FOREWORD 


i  This  standard  defines  the  Government’s  requirements  and  expectations  for  contractor 

peiformance  in  defense  system  acquisitions  and  technology  developments. 

This  new-issue  SMC  standard  comprises  the  text  of  The  Aerospace  Corporation  report 
number  TOR-2007(8506)-7164,  entitled  Requirements  for  End-Of-Life  Disposal  of  Satellites 
Operating  in  Orbits  with  Perigees  below  2000  Kilometers. 

3 .  Beneficial  comments  (recommendations,  changes,  additions,  deletions,  etc.)  and  any 
pertinent  data  that  may  be  of  use  in  improving  this  standard  should  be  forwarded  to  the  following 
addressee  using  the  Standardization  Document  Improvement  Proposal  appearing  at  the  end  of  this 
document  or  by  letter: 


Division  Chief,  SMC/EAE 
SPACE  AND  MISSILE  SYSTEMS  CENTER 
Air  Force  Space  Command 
483  N.  Aviation  Blvd. 

El  Segundo,  CA  90245 

4.  This  standard  has  been  approved  for  use  on  all  Space  and  Missile  Systems 
Center/Air  Force  Program  Executive  Office  -  Space  development,  acquisition,  and 
sustainment  contracts. 


DAVID  E.  SWANSON,  COL,  USAF 
SMC  Chief  Engineer 
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1.  Scope 


According  to  the  U.S.  Government  Orbital  Debris  Standard  Practices, 1  a  spacecraft  or  upper  stage 
with  perigee  altitude  below  2000  km  in  its  final  mission  orbit  may  be  disposed  of  by  using  an 
atmospheric  reentry  option.  The  guideline  for  atmospheric  reentry  is  stated  as  follows: 

“Leave  the  structure  in  an  orbit  in  which,  using  conser\’ative  projections  for  solar  activity, 
atmospheric  drag  will  limit  the  lifetime  to  no  longer  than  25  years  after  completion  of  mission.  If 
drag  enhancement  devices  are  to  be  used  to  reduce  the  orbit  lifetime,  it  should  be  demonstrated 
that  such  devices  will  significantly  reduce  the  area-time  product  of  the  system  or  will  not  cause 
spacecraft  or  large  debris  to  fragment  if  a  collision  occurs  while  the  system  is  decaying  from 
orbit.  ’’ 

The  Inter-Agency  Space  Debris  Coordination  Committee  (IADC)  Space  Debris  Mitigation 
Guidelines2  state: 

“Whenever  possible  space  systems  that  are  terminating  their  operational  phases  in  orbits  that  pass 
through  the  LEO  region,  or  have  the  potential  to  interfere  with  the  LEO  region,  should  be  de-orbited 
(direct  re-entry  is  preferred)  or  where  appropriate  maneuvered  into  an  orbit  with  a  reduced  lifetime. 
Retrieval  is  also  a  disposal  option. 

A  space  system  should  be  left  in  an  orbit  in  which,  using  an  accepted  nominal  projection  for  solar 
activity,  atmospheric  drag  will  limit  the  orbital  lifetime  after  completion  of  operations.  A  study  on  the 
effect  of  post-mission  orbital  lifetime  limitation  on  collision  rate  and  debris  population  growth  has 
been  performed  by  the  IADC.  This  IADC  and  some  other  studies  and  a  number  of  existing  national 
guidelines  have  found  25  years  to  be  a  reasonable  and  appropriate  lifetime  limit.  If  a  space  system  is 
to  be  disposed  of  by  re-entry  into  the  atmosphere,  debris  that  survives  to  reach  the  surface  of  the 
Earth  should  not  pose  an  undue  risk  to  people  or  property.  This  may  be  accomplished  by  limiting  the 
amount  of  surviving  debris  or  confining  the  debris  to  uninhabited  regions,  such  as  broad  ocean 
areas.  Also,  ground  environmental  pollution,  caused  by  radioactive  substances,  toxic  substances  or 
any  other  environmental  pollutants  resulting  from  on-board  articles,  should  be  prevented  or 
minimized  in  order  to  be  accepted  as  permissible. 

In  the  case  of  a  controlled  re-entry  of  a  space  system,  the  operator  of  the  system  should  inform  the 
relevant  air  traffic  and  maritime  traffic  authorities  of  the  re-entry  time  and  trajectory >  and  the 
associated  ground  area.  ” 

This  report  specifies  requirements  for  (1)  planning  for  disposal  of  satellites  operating  in  low  Earth 
orbits  with  perigee  altitude  below  2000  km,  (2)  selecting  the  final  disposal  orbits  whose  perigee 
altitude  is  low  enough  to  ensure  a  natural  decay  due  to  atmospheric  drag  within  25  years,  (3) 
executing  the  disposal  maneuver  successfully,  and  (4)  depleting  all  energy  sources  onboard  the 
vehicle  prior  to  end-of-life  to  minimize  the  possibility  of  a  debris-producing  event.  Techniques  for 
planning  and  determining  the  final  disposal  orbit  are  provided,  considering  current  internationally 
accepted  guidelines  and  operational  best  practices. 

2.  Application 

These  requirements  shall  be  applied  to  the  end-of-life  disposal  of  satellites  operating  in  LEO  (Low 
Earth  Orbit)  or  HEO  (Highly  Elliptical  Orbit)  with  low  perigee  altitudes  (Hp  <  2000  km).  The 
requirements  are  applicable  to  all  DOD/NASA/NOAA/civil  space  programs  with  perigee  altitude 
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lower  than  2000  km,  such  as  DMSP,  NPOESS,  Iridium,  Molniya,  and  GTO  (Geosynchronous 
Transfer  Orbit)  rocket  bodies. 

3.  Terms  and  Definitions 


•  End  of  Life  (EOL):  The  point  in  time  at  which  a  satellite  is  permanently  turned  off 
(nominally  after  the  space  system  has  been  placed  in  a  disposal  orbit)  or  when  control  of  the 
satellite  from  the  ground  is  no  longer  possible. 

•  End  of  Mission:  End  of  mission  refers  to  the  completion  of  the  scheduled  mission  of  the 
space  system,  or  when  the  mission  is  halted  either  because  of  a  failure  in  space  or  a  voluntary 
decision. 

•  LEO  Region:  The  LEO  region  is  a  spherical  region  that  extends  from  the  Earth’s  surface  up 
to  an  altitude  of  2,000  km.  . 

•  Drag  Enhancement  Device:  A  drag  enhancement  device  on  a  spacecraft  effectively 
increases  the  cross-sectional  area  of  the  spacecraft  and  thereby  increases  the  atmospheric  drag 
effect  to  shorten  the  decay  time.  One  such  effective  device  is  the  balloon. 

•  Direct  Reentry:  Without  using  the  atmospheric  drag  decay,  the  spacecraft  is  inserted  into  an 
orbit  with  perigee  low  enough  to  ensure  a  ground  impact  after  a  half  revolution  of  the  orbit. 


4.  Symbols  and  Acronyms 


a 

A 

* p 

P 

CD 

e 

EOL 
FI  0.7 


F  10.7  Bar 
GEO 
Ha 
HP 

IADC 

ISO 

LEO 

m 

RAAN 


Re 


Orbit  semimajor  axis 

Satellite  cross-sectional  area  with  respect  to  the  relative  wind 

Earth  daily  geomagnetic  index 

Ballistic  coefficient  of  satellite  =  CD  *  A  / m 

Satellite  drag  coefficient 

Orbit  eccentricity 

End  of  life 

Solar  radio  flux  observed  daily  at  2800  MHz  (10.7  cm)  in  solar  flux  units 
(10'22W  m'2  Hz'1) 

Solar  radio  flux  at  2800  MHz  (10.7  cm),  averaged  over  three  solar  rotations 

Geosynchronous  Earth  orbit 

Apogee  altitude  =  a  (1  +  e)  -  Rg 

Perigee  altitude  =  a  (1  -  e)  -  Rg 

Inter-Agency  Space  Debris  Coordination  Committee 

International  Standard  Organization 

Low  Earth  orbit 

Mass  of  satellite 

Orbit  Right  Ascension  of  the  Ascending  Node;  the  angle  between  the  vernal  equinox 
and  the  orbit  ascending  node,  measured  counter-clockwise  in  the  equatorial  plane, 
looking  in  the  -Z  direction. 

Equatorial  radius  of  the  Earth  =  6378. 137  km 
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5.  Primary  Requirements 

5.1  Disposal  Maneuver  Planning 

An  “End  of  Mission  Disposal  Plan”  (EOMDP)  shall  be  developed,  maintained,  and  updated  in  all 
phases  of  mission  and  spacecraft  design  and  operations.  An  accurate  estimate  of  fuel  shall  be  required 
to  ensure  orbit  decay  within  25  years.  This  document  shall  include: 

•  Details  of  the  nominal  mission  orbit 

•  Details  of  the  targeted  disposal  orbit 

•  Estimates  of  the  fuel  required  for  the  disposal  action 

•  Identity  of  system  and  capabilities  required  for  successful  completion  of  the  disposal  action 

•  Criteria  that  dictate  initiation  of  the  disposal  action  when  met 

•  Identities  of  energy  sources  that  must  be  depleted  prior  to  end-of-life 

•  A  timeline  for  initiation  and  execution  of  the  disposal  action 

•  A  timeline  for  depletion  of  remaining  energy  sources 

•  A  list  of  those  individuals  and/or  entities  to  be  notified  of  the  end-of-mission  and  disposal  and 
a  timeline  for  notification 

5.2  Reliability  for  Disposal 

The  space  system  shall  be  designed  such  that  the  probability  of  successful  end-of-mission  disposal, 
including  depletion  of  energy  sources,  equals  or  exceeds  0.9  at  end  of  design  mission  life.  The 
disposal  success  probability  shall  be  evaluated  as  conditional  probability  (weighted  on  the  mission 
success).  Appendix  E  provides  a  discussion  of  the  conditional  probability.  Details  of  the  design  that 
provide  the  basis  for  the  probability  estimate  shall  be  included  in  the  EOMDP. 

5.3  Criteria  for  Executing  Disposal  Action 

Specific  criteria  for  initiation  of  the  disposal  action  shall  be  developed,  included  in  the  EOMDP,  and 
monitored  throughout  the  mission  life.  Projections  of  mission  life  based  on  these  criteria  shall  be  a 
regular  part  of  mission  status  reviews  during  the  mission.  The  status  of  these  criteria  shall  be 
presented  at  periodic  (at  least  yearly)  mission  reviews. 

5.4  Energy  Depletion  Planning 

Independent  of  the  success  or  failure  of  other  aspects  of  the  disposal  action,  a  plan  shall  be  developed 
to  deplete  all  energy  sources,  and  safe  the  vehicle  in  accordance  with  ISO  passivation  standards 
before  the  final  demise  of  the  spacecraft.  (Passivation  or  safing  is  the  elimination  of  all  stored  energy 
on  a  spacecraft  or  orbital  stages  to  reduce  the  chance  of  break-up.  Typical  passivation  measures 
include  venting  or  burning  excess  propellant,  discharging  batteries,  and  relieving  pressure  vessels.) 
The  objective  shall  be  to  ensure  that  actions  necessary  to  safe  the  vehicle  are  taken  before  critical 
systems  are  lost.  The  plan  shall  include  criteria  that  define  when  the  safing  actions  are  to  be  taken,  the 
rationale  for  each  criterion,  and  a  schedule  for  safing  actions.  The  plan  shall  be  included  in  the 
EOMDP. 
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6.  Disposal  Planning  Requirements 


Planning  activities  for  end-of-mission  disposal  shall  start  in  mission  design.  Planning  for  the  actual 
disposal  action  shall  begin  at  least  six  months  prior  to  the  date  of  re-orbit  maneuvers.  The  following 
steps  shall  be  followed  in  all  mission  phases  and  shall  be  documented  in  the  EOMDP. 

6.1  Estimate  Ballistic  Coefficient  (CdA/m)  Considering  Tumbling 

The  first  step  in  planning  for  the  disposal  of  a  LEO  or  EIEO  spacecraft  is  to  estimate  the  ballistic 
coefficient  (CdA/m)  considering  spacecraft  tumbling.  It  is  generally  believed  that  the  average  value 
between  the  maximum  and  minimum  cross-section  area  should  be  used  to  compute  the  ballistic 
coefficient  of  a  tumbling  spacecraft.  The  proper  value  of  the  dimensionless  drag  coefficient,  Cd,  is  2 
for  a  typical  spacecraft.  A  discussion  on  the  drag  coefficient  is  given  in  Appendix  A. 

6.2  Determine  Orbit  Lifetime  of  Nominal  Mission  Orbit 

The  second  step  is  to  predict  the  orbit  lifetime  of  the  nominal  mission  orbit  beginning  at  the  end-of- 
life.  No  further  planning  is  required  if  the  predicted  orbit  lifetime  based  on  one  of  the  following 
recommended  methods  is  within  25  years.  See  Section  6.3  if  the  predicted  lifetime  is  longer  than  25 
years. 

To  predict  orbit  lifetime: 

•  Perform  long-term  orbit  propagations  using  either  semi-analytical  or  numerical  integration 
tools.  Choose  the  up-to-date  atmospheric  density  model,  such  as  JB2006  or  MSISEOO.  The 
recommended  long-term  propagation  tools  and  the  solar  flux  data  are  discussed  in  detail  in 
Section  7  and  8. 

•  Determine  orbit  lifetime  from  the  empirical  equations  provided  in  Appendix  B.  These 
empirical  equations  are  not  applicable  to  the  three  special  orbit  types:  sun-synchronous  orbits, 
Molniya  orbits  and  Magic  orbits.  These  orbits  are  discussed  in  Appendix  D.  The  long-term 
lifetime  prediction  should  be  performed  using  semi- analytical  or  numerical  integration  tools 
as  in  a)  above. 

6.3  Determine  the  Strategy  for  Ensuring  25-Year  Lifetime 

There  are  two  basic  strategies  to  ensure  a  25-year  lifetime  if  the  nominal  mission  orbit  lifetime  is 
longer  than  25  years.  The  first  strategy  is  to  perform  maneuvers  to  lower  the  orbit  altitude  to  achieve 
the  25-year  lifetime.  The  second  strategy  is  to  use  drag  enhancement  devices,  such  as  balloons  or 
parachutes. 

6.3.1  Lower  Orbit  Altitude 

The  recommended  strategy  to  lower  the  orbit  altitude  is  to  perform  a  single-bum  at  apogee  to  lower 
the  perigee  altitude  to  the  predetermined  value  that  ensures  a  25-year  lifetime.  It  has  been  proven 
through  numerical  simulations  that  a  single-bum  maneuver  requires  less  propellant  than  a  two-bum 
Hohmann  type  transfer  to  ensure  a  lifetime  less  than  25  years. 

The  procedure  to  determine  the  amount  of  AV  required  is  through  either  numerical  orbit  propagation 
or  table  lookup,  as  discussed  in  section  6.2  above.  The  single  bum  should  be  executed  near  a 
predetermined  time  such  that  the  new  argument  of  perigee  is  at  or  close  to  270  deg.  This  optimum 
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eccentricity  vector3  can  shorten  the  25-year  lifetime  by  1  to  20  percent,  depending  on  initial 
eccentricity  and  perigee  altitude.  Consequently,  the  amount  of  required  AV  for  ensuring  decay  within 
25  years  may  be  significantly  reduced. 

6.3.2  Use  Drag  Enhancement  Devices 

This  strategy  requires  the  use  of  drag  enhancement  devices,  such  as  balloons  or  parachutes.  The 
average  ballistic  coefficient  of  the  spacecraft  with  drag  enhancement  device  deployed  must  be 
estimated.  The  new  ballistic  coefficient  (CDA/m)  is  determined  by  iteration  using  either  one  of  the 
two  methods  described  in  Section  6.2.  This  strategy  may  also  require  an  apogee  bum  to  lower  the 
perigee,  depending  on  the  effectiveness  of  the  drag  device.  In  this  strategy,  the  deployment  of  the 
drag  device  and  the  apogee  bum,  if  needed,  should  also  be  timed  as  discussed  in  a)  for  achieving  the 
initial  optimum  eccentricity  vector.  An  assessment  of  the  drag  enhancement  device  is  recommended 
to  demonstrate  that  either  the  device  will  significantly  reduce  the  collision  risk  of  the  system  by  the 
reduction  of  the  orbit  lifetime  or  that  the  device  will  not  cause  spacecraft  or  large  debris  to  fragment  if 
a  collision  occurs  while  the  system  is  decaying  from  orbit. 

6.4  Develop  Basic  Maneuver  Requirements  for  Achieving  Disposal  Orbit 

The  mission  operations  shall  develop  maneuver  sequences  to  achieve  the  predetermined  disposal  orbit 
and  to  assure  that  drag  enhancement  devices  are  properly  deployed  if  this  strategy  is  used.  The 
predetermined  disposal  orbit  shall  satisfy  the  following  conditions: 

•  The  perigee  altitude  must  be  low  enough  to  ensure  that  the  spacecraft  shall  decay  through 
atmospheric  drag  within  25  years. 

•  The  argument  of  perigee  shall  be  either  close  to  270  deg  or  the  estimated  value  of  the 
optimum  eccentricity  vector. 

In  all  cases,  the  space  system  shall  be  passivated  prior  to  end  of  life. 

7.  Orbit  Lifetime  Assessment 

The  orbit  lifetime  estimation  process  is  represented  generically  in  Figure  1. 


Spacecraft  Orbit  Initial 
Conditions  at  Epoch 


Spacecraft  Ballistic 
Characteristics,  Attitude 
Rules 


Anticipated  Long- 
Duration  Solar  & 
Geomagnetic  Activity 


Orbit  Integration  and 
Propagation  Package 


Orbit  Lifetime 
Estimate 


Atmosphere  Model 


Figure  1.  Orbit  lifetime  assessment  process. 
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7.1  Orbit  Lifetime  Assessment  Methods  and  Applicability 

There  are  three  basic  analysis  methods  used  to  estimate  orbit  lifetime.  Method  1  is  simply  a  table 
lookup,  graphical  analysis,  or  evaluation  of  equations  fit  to  pre-computed  orbit  lifetime  estimation 
data  obtained  via  the  extensive  and  repetitive  application  of  Methods  2  and/or  3.  Method  2  utilizes  a 
definition  of  mean  orbital  elements4,5  semi-analytic  orbit  theory  and  average  satellite  ballistic 
coefficient  to  permit  the  very  rapid  integration  of  the  equations  of  motion,  while  still  retaining 
reasonable  accuracy.  Method  3,  clearly  the  highest  fidelity  model,  utilizes  a  numerical  integrator  with 
a  detailed  gravity  model,  third-body  effects,  solar  radiation  pressure,  and  a  detailed  satellite  ballistic 
coefficient  model. 

7.1.1  Method  1:  Numerical  Table  Look-up,  Analysis  and  Fit  Equation  Evaluations 

In  this  first  method,  one  uses  tables,  graphs  and  equations  representing  data  that  was  generated  by 
exhaustively  using  Methods  2  and  3  (see  below).  The  tables,  graphs  and  equations  provided  in  this 
report  can  help  the  analyst  assess  orbit  lifetime  for  their  particular  case  of  interest;  all  such  Method  1 
data  in  this  report  were  generated  using  Method  2  approaches. 

7.1.2  Method  2:  Rapid  Semianalytical  Orbit  Propagation 

Method  2  analysis  tools  utilize  semianalytic  propagation  of  mean  orbit  elements4, 5  influenced  by 
gravity  zonals  J2  and  J3  and  selected  atmosphere  models.  The  primary  advantage  of  this  approach 
over  direct  three-dimensional  numerical  integration  of  the  equations  of  motion  (Method  3)  is  that 
long-duration  orbit  lifetime  cases  can  be  quickly  analyzed  (e.g.,  5  seconds  CPU  time  for  a  30-year 
orbit  lifetime  case).  While  incoiporation  of  an  attitude-dependent  ballistic  coefficient  is  doable  for 
this  method,  an  average  ballistic  coefficient  is  typically  used. 

The  Method  2  results  provided  in  this  report  were  generated  using  the  LIFETIME  program,6, 7’ 8 
developed  and  maintained  by  The  Aerospace  Corporation,  and  the  QPROP  program,  developed  and 
maintained  by  1  Earth  Research,  LLC.  Both  tools  are  PC-based,  use  Method  2-type  orbit  lifetime 
estimators,  using  semi-analytic  propagation  of  mean  orbit  elements  coupled  with  gravity  zonals  J2  and 
J3  and  selected  atmosphere  models  (including  MSISE2000,  Jacchia-Bowman  2006,  Jacchia  1971, 
etc.).  The  LIFETIME  analysis  tool  has  been  used  to  analyze  orbit  lifetime  and  satellite  reentry  by 
various  Government  and  industrial  organizations.  Its  accuracy  has  been  validated  by  high-precision 
numerical  integration  results  (Method  3  type)  as  well  as  actual  satellite  decay  data. 

The  LIFETIME  and  QPROP  software  tools  were  exercised  in  this  study  to  jointly  analyze  over  eight 
million  orbit  lifetime  cases,  spanning  a  variety  of  times-into-the-solar-cycle,  inclinations,  perigee 
altitudes  (Hp),  apogee  altitudes  (Ha),  and  ballistic  coefficients. 

7.1.3  Method  3:  High-precision  numerical  integration 

Method  3  is  the  direct  numerical  integration  of  all  accelerations  in  Cartesian  space,  with  the  ability  to 
incorporate  a  detailed  gravity  model  (e.g.,  using  a  larger  spherical  harmonics  model  to  address 
resonance  effects),  third-body  effects,  solar  radiation  pressure,  vehicle  attitude  rules  or  aero-torque- 
driven  attitude  torques,  and  a  detailed  satellite  ballistic  coefficient  model  based  on  the  angle-of-attack 
with  respect  to  the  relative  wind.  Atmospheric  rotation  at  the  Earth’s  rotational  rate  is  also  easily 
incorporated  in  this  approach.  The  only  negative  aspects  to  such  simulations  are  (1)  they  run  much 
slower  than  Method  2;  (2)  many  of  the  detailed  data  inputs  required  to  make  Method  3  realize  its  full 
accuracy  potential  are  not  available;  and  (3)  any  gains  in  orbit  lifetime  prediction  accuracy  are 
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frequently  overwhelmed  by  inherent  inaccuracies  of  atmospheric  modeling.  However,  this  is 
undoubtedly  the  most  accurate  method  to  analyze  a  few  select  cases  where  such  detailed  model  inputs 
are  known. 

8.  Atmospheric  Density  Modeling 

The  three  biggest  factors  in  orbit  lifetime  estimation  are  (1)  the  selection  of  an  appropriate 
atmosphere  model  to  incorporate  into  the  orbit  acceleration  formulation;  (2)  the  selection  of 
appropriate  inputs  to  such  a  model;  and  (3)  ballistic  coefficient  determination. 

8.1  Atmospheric  Drag  Models 

A  wide  variety  of  atmosphere  models  are  available  to  the  orbit  analyst.  The  background,  technical 
basis,  utility  and  functionality  of  these  atmosphere  models  are  described  in  detail  in  Ref.  9.  This 
report  does  not  presume  to  dictate  which  atmosphere  model  the  analyst  shall  use.  However,  it  is 
worth  noting  that  in  general,  the  heritage,  expertise  and  especially  the  observational  data  that  went 
into  creating  each  atmosphere  model  play  a  key  role  in  that  model’s  ability  to  predict  atmospheric 
density,  which  is  in  turn  a  key  factor  in  estimating  orbit  lifetime.  Many  of  the  early  atmosphere 
models  were  low  fidelity  and  were  created  on  the  basis  of  only  one,  or  perhaps  even  just  a  part  of  one, 
solar  cycle’s  worth  of  data. 

The  advantage  of  some  of  these  early  models  is  that  they  typically  run  much  faster  than  the  latest 
high-fidelity  models  (see  Figure  2)  without  losing  a  great  deal  of  accuracy.  However,  the  use  of 
atmosphere  models  that  were  designed  to  fit  a  select  altitude  range  (e.g.,  the  “exponential” 
atmosphere  model  shown  in  Figure  2)  or  models  that  do  not  accommodate  solar  activity  variations 
should  be  avoided,  as  they  miss  too  many  atmospheric  density  variations  to  be  sufficiently  accurate. 

There  are  some  early  models  (e.g.,  the  Jacchia  1 97 1 10)  that  do  accommodate  solar  activity  variations 
yet  run  very  fast;  these  models,  while  perhaps  not  ideal,  can  work  well  for  long-duration  orbit  lifetime 
studies  where  numerous  cases  are  to  be  examined.  A  crude  comparison  of  a  sampling  of  atmosphere 
models  for  a  single  test  case  is  shown  in  Figures  3  and  4,  illustrating  the  range  of  temperatures  and 
densities  exhibited  by  the  various  models. 

Use  of  the  more  recent  atmosphere  models  (e.g.,  MS1SE200011  and  Jacchia-Bowman12)  are 
encouraged  because  they  have  substantially  more  atmospheric  drag  data  incorporated  as  the 
foundation  of  their  underlying  assumptions.  A  recent  study  by  Marcos,  et  al.13,  indicates  that  the 
standard  deviations  for  the  JB2006  model  are  lower  at  all  altitudes.  Both  the  improved  solar  indices 
and  semiannual  formulations  contribute,  approximately  equally,  to  the  successful  increase  in  JB2006 
model  precision. 
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Figure  2.  Atmosphere  model  runtime  (1M  evals). 
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Figure  3.  Temperature  comparison  by  atmosphere  model. 
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Figure  4.  Comparison  of  a  small  sampling  of  atmosphere  models. 

8.2  Long-Term  Solar  Flux  and  Geomagnetic  Index  Prediction 

The  goal  of  the  standard  is  to  help  minimize  space  object  overpopulation  for  objects  in  the  LEO- 
crossing  orbit  regimes.  This  is  accomplished  by  standardizing  on  an  anticipated  post-mission  lifetime 
of  25  years  or  less.  Utilization  of  the  higher- fidelity  atmosphere  models  mentioned  in  the  previous 
section  requires  the  orbit  analyst  to  specify  the  solar  and  geomagnetic  indices  required  by  such 
models.  As  can  be  seen  by  the  careful  comparison  of  Figures  5  through  8,  attempts  to  ‘generalize’  or 
abstract  solar  indices  are  fraught  with  error,  since  they  do  not  capture  the  true  upper  and  lower 
variations.  Also,  note  that  the  “estimated”  values  from  typical  solar  cycle  predictions  appear  to  be 
located  at  the  mean  value  of  the  upper  and  lower  curves.  But  as  shown  in  Figure  6,  the  frequency  of 
occurrence  is  highest  near  the  lowest  prediction  boundary. 


There  are  other  problems  with  the  use  of  the  F10.7  curves  shown  in  Figure  5,  notably  that  the  timing 
of  the  cycle  minima  and  maxima  are  always  imprecise.  The  resultant  time  bias  that  such  a  prediction 
error  would  introduce  could  easily  yield  large  FI 0.7  prediction  errors  of  100%  or  more. 

The  long  time  duration  currently  being  advocated  by  the  IADC  (25  years,  although  this  is  still  under 
review)  suggests  a  different  approach  for  the  standard.  Rather  than  attempt  to  predict  a  ‘mean’  curve, 
which  would  ignore  the  highly  nonlinear  aspects  of  solar  storms  and  quiet  periods,  the  adopted 
approach  is  to  utilize  historical  F10.7  and  Ap  (geomagnetic)  data14' 15  directly,  thereby  retaining  the 
wide  variability  and  nonlinearity  of  the  data.  Care  is  taken  to  adopt  the  adjusted  F10.7  values  (adjusted 
to  correct  for  Earth-sun  distance  variations).  Note  that  we  already  have  more  than  five  solar  cycles  of 
data  to  choose  from.  Processing  of  this  data  maps  each  coupled  and  correlated  triad  of  datum  (F10.7, 
F10.7  Bar,  and^4Q  into  a  single  solar  cycle  range  of  10.82546  years  (3954  days). 

Implementation  of  this  historical  data,  mapped  into  a  single  solar  cycle  (Figures  9  through  11),  then 
uses  random  draws  upon  the  number  of  data  triads  available  on  any  day  within  the  single  solar  cycle. 
Since  we  have  accumulated  daily  data  since  14  February  1947,  on  any  given  day  within  the  3954-day 
solar  cycle  we  have  at  least  five  data  triads  to  choose  from.  It  is  important  that  the  random  draw 
retain  the  integrity  of  each  data  triad,  since  F10.7,  F10.7  Bar,  and  Ap  are  interrelated. 

In  summary,  the  selected  method  used  in  the  QPROP  analysis  tool  for  modeling  solar  and 
geomagnetic  data  is  to  select  a  new  coupled  triad  of  (F10.7,  F10.7  Bar,  and  Ap)  data  for  each  day  of  the 
orbit  lifetime  simulation,  thereby  simulating  the  drag  effect  resulting  from  solar  and  geomagnetic 
variations  consistent  with  historical  trends  for  these  data.  The  C++  code  that  QPROP  uses  to 
implement  this  atmospheric  variation  strategy  has  been  made  publicly  available  so  that  users  of  the 
proposed  orbit  lifetime  standard  can  adopt  a  standardized  F10.7,  F10.7  Bar,  and^,  prediction  if  desired. 
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Figure  5.  Solar  flux  estimated  upper,  lower  and  mean. 


Figure  6.  Distribution  of  F10.7  deviates. 
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Figure  7.  F10.7  Solar  data. 
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Figure  8.  Averaged  FI 0.7  Solar  data  with  minima  &  maxima. 
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FI 0.7  Bar  Solar  Flux  Index 
(Avg  F10.7  over  3  Solar  Rotations  ~  81  days) 


F10.7  Normalized  to  Avg  Cycle 


Figure  9.  F10.7  Normalized  to  average  cycle. 
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Figure  10.  F10.7  Bar  normalized  to  average  cycle. 
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Ap  Normalized  to  Avg  Cycle 


Days  From  Start  of  Cycle 


Figure  1 1.  Ap  normalized  to  average  cycle. 


It  can  also  be  seen  from  Figure  1 1  that^,  is  (1)  unpredictable;  (2)  loosely  correlated  with  the  solar 
cycle;  and  (3)  volatile.  Figure  12  demonstrates  that  density  varies  greatly  depending  uponJp;  thus,  a 
geomagnetic  storm  could  induce  large  decreases  in  orbital  energy  (orbit  decay)  during  such  storms 
that  the  use  of  some  average  Ap  value  would  miss.  But  while  this  dependency  does  exist,  the 
unpredictable  nature  of  Ap  shown  in  Figure  1 1  again  illustrates  that  perhaps  the  best  implementation 
of  geomagnetic  index  modeling  is  the  use  of  random  draws  of  historical  data,  such  as  we  adopt  in  this 
report  for  orbit  lifetime  standards  purposes. 
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Figure  12.  Log  (density)  variation  as  a  function  of^  value. 
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9.  Conclusions 


The  process  of  orbit  lifetime  estimation  has  been  examined,  with  the  goal  of  standardizing  the  space 
community’s  methods  for  handling  such  uncertain  parameters  as  solar  and  geomagnetic  indices, 
ballistic  coefficient  prediction,  and  thermospheric  cooling.  The  two  primary  methods  of  assessing 
orbit  lifetime  (Methods  2  and  3),  together  with  their  derived  analysis  products  (Method  1),  provides 
spacecraft  designers,  owners,  and  operators  with  a  variety  of  assessment  methods  that  span  a  wide 
range  of  orbit  initial  conditions,  implementation  effort  required,  and  utility. 

Pertaining  to  the  IADC  original  25 -year  guideline  for  post- mission  residency  in  the  0-2000  km  LEO 
band,  numerous  computer  runs  have  been  accomplished  that  show  the  family  of  post-safmg  orbit 
initial  conditions  that  were  found  to  exhibit  a  25-year  post-mission  lifetime.  For  a  user-provided 
perigee  altitude,  it  is  logical  that  if  the  user’s  apogee  altitude  is  below  this  25-year  proximity  region, 
then  the  user’s  post-mission  orbit  will  decay  in  less  than  25  years  and  will  meet  IADC  guidelines  as 
well  as  the  proposed  ISO  standard.  Conversely,  if  for  the  same  user-provided  perigee  altitude  the 
apogee  altitude  is  above  the  25-year  proximity  region,  then  the  user’s  orbit  will  not  decay  within  25 
years  and  therefore  violates  the  IADC  guideline  and  ISO  standard.  Finally,  if  the  user’s  specified 
orbit  falls  within  the  region  of  the  25-year  constraint,  then  it  may  or  may  not  be  an  issue;  perhaps 
using  increased  modeling  fidelity  (e.g.,  switching  to  Method  3,  coupled  with  the  statistical  solar  and 
geomagnetic  indices)  would  allow  the  analyst  to  make  a  better  determination  as  to  whether  the  IADC 
guideline  and  ISO  standard  requirements  will  be  met. 
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Appendix  A:  Estimating  Ballistic  Coefficient  (CDA/m)  Considering  Tumbling 


The  first  step  in  planning  for  a  LEO  or  HEO  spacecraft  disposal  is  to  estimate  the  ballistic  coefficient 
(CDA/m).  Cross-sectional  area  and  drag  coefficient  estimations  are  examined  separately.  The 
spacecraft  mass  is  typically  assumed  to  be  constant  from  EOL  until  orbit  decay. 

A.1  Estimating  Cross-Sectional  Area  with  Tumbling  and  Stabilization 

The  average  cross-sectional  area  should  be  used  to  compute  the  ballistic  coefficient  of  a  tumbling 
spacecraft.  This  can  be  as  simple  as  a  straightforward  two-point  average  of  the  minimum  and 
maximum  cross-sectional  areas,  or  this  approach  can  be  extensively  refined  by  integrating  the  cross- 
sectional  area  of  the  spacecraft  across  all  anticipated  tumbling  attitudes,  and  then  dividing  the  result 
by  the  difference  between  the  limits  of  integration.  The  analyst  is  then  left  with  a  properly  weighted 
average  cross-sectional  area.  For  satellites  with  a  large  length-to-diameter  ratio,  the  analyst  should 
consider  whether  gravity  gradient  stabilization  will  occur  and  adjust  the  cross-sectional  area 
accordingly.  Similarly,  for  satellites  that  have  a  large  aero-torque  moment  (i.e.,  the  center-of-gravity 
and  center-of-pressure  are  suitably  far  apart  and  the  aerodynamic  force  is  suitably  large),  the  analyst 
should  consider  whether  the  satellite  would  experience  drag-induced  passive  attitude  stabilization  and 
adjust  the  cross-sectional  area  accordingly. 

A.2  Estimating  Drag  Coefficient 

The  proper  value  of  the  dimensionless  drag  coefficient,  CD,  is  2  for  a  typical  spacecraft. 

The  atmospheric  drag  force  is  directly  proportional  to  the  ballistic  coefficient,  CDA/m,  contained  in 
the  equations  of  motion.  The  drag  coefficient,  CD,  depends  on  the  shape  of  the  satellite  and  the  way 
air  molecules  collide  with  it.  In  free  molecular  flow,  CD  is  usually  assumed  to  be  the  theoretical  value 
of  2.  However,  for  certain  geometric  configurations  such  as  spheres,  cylinders  and  cones,  the  value  of 
axial  drag  coefficient,  CD,  can  be  evaluated  more  precisely  than  previously  noted  provided  something 
is  known  about  the  flow  regime  and  reference  area. 

As  discussed  in  Ref.  16,  the  atmosphere  above  150  km  can  be  represented  by  free  molecular  flow 
conditions.  The  free  molecular  flow  regime  is  defined  by  the  Knudsen  number,  Kn,  which  is  the  ratio 
of  the  mean  free  path  of  a  molecule  to  some  reference  length  of  the  body  pertinent  to  the  flow  field. 
For  example,  the  mean  free  path  at  150  km  altitude  is  about  6.5  m.  Therefore,  for  a  1.6-m  diameter 
spacecraft,  the  Knudson  number  is  about  4,  which  is  in  the  transition-to-free-molecular  flow  region. 
The  drag  coefficient  of  simple  symmetrical  bodies,  such  as  spheres,  cylinders  and  cones,  may  be 
analytically  expressed18  as  a  function  of  spacecraft  to  molecule  speed  ratio  s: 

s  =  V  /  vmp  (A- 1 ) 

where  V  is  the  spacecraft’s  velocity  and  vmp  is  the  most  probable  molecular  speed.  The  mean 
molecular  speed,  v*,  is  given  in  most  atmospheric  tables  such  as  ARDC  1959  and  can  be  corrected  to 
given  vmp  as  follows: 


vmp  =  0.8862  v*  (A-2) 

For  example,  assume  a  satellite  is  at  555  km  with  a  velocity  of  about  7.6  km/sec.  The  most  probable 
molecular  speed  is  about  1.19  km/sec.  Thus,  the  molecular  speed  ratio  s  is  about  6.4.  At  740  km 
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altitude,  s  =  5.9.  Figure  A-l  gives  the  drag  coefficients  for  a  sphere  for  both  diffuse  and  specular 
reflection16  (i.e.,  angle  of  incidence  equals  angle  of  reflection).  For  design  or  analysis  purposes, 
curves  for  diffuse  reflection  should  be  used,  since  the  fraction  of  diffusely -reflected  particles  is  quite 
high  for  free  molecular  flow.  Figure  A-2  gives  drag  coefficients  for  a  cylinder16. 

Note  that  for  high  molecular  speed  ratios,  Cd  =  2  is  reached  for  spherical  bodies  and  CD  =  3  for 
cylindrical  bodies  with  the  relative  wind  along  the  axis  of  rotation  of  the  cylinder.  These  values  are 
quite  adequate  as  a  first  approximation.  For  speed  ratios  much  lower  than  5  for  spherical  bodies  or 
much  lower  than  2  for  cylindrical  bodies,  the  above  approximation  is  no  longer  good,  as  shown  in 
Figures  A-l  and  A-2.  For  example,  in  the  LIFETIME  program  the  CD  value  is  computed  from  an 
empirical  equation  as  a  function  of  altitude  determined  from  Newtonian  flow  when  the  altitude  is 
lower  than  125  km.  A  theoretical  discussion17  of  drag  coefficients  in  1965  also  supports  the  above 
approximation  (CD  =  2.2  for  spheres). 


Specular 

Diffuse 


Molecular  speed  ratio,  s 

Figure  A-l.  Sphere  drag  coefficient,  diffuse  and  specular  reflection. 


Specula 

Diffuse 


Figure  A-2.  Cylinder  drag  coefficient,  diffuse  and  specular  reflection. 
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Based  on  the  approximate  mass,  shape,  size  and  orientation  of  various  LEO  satellites,  values  of 
ballistic  coefficients  and  drag  coefficients  were  estimated  from  post-flight  orbit  determination.18  The 
typical  upper  and  lower  limits  of  drag  coefficients  are  4  and  2,  respectively.18  Table  A-l  shows  the 
maximum  and  minimum  values  of  the  16  satellites  in  low  Earth  orbit.  The  estimated  drag  coefficients 
are  bounded  between  2  and  4. 

Table  A-l.  Estimated  Ranges  of  Drag  Coefficients  for  Low-Earth  Orbit  Satellites* 


Satellite 

Mass,  kg 

Shape 

Max  XA 
(m2) 

Min  XA 
(m2) 

Max 

XA 

CD 

Min  XA 

cD 

Type  of 
mission 

Oscar-1 

5 

box 

0.075 

0.0584 

4 

2 

Comm. 

lntercos-16 

550 

cylinder 

2.7 

3.16 

2.67 

2.1 

Scientific 

Viking 

277 

octagon 

2.25 

0.833 

4 

2.6 

Scientific 

Explorer-1 1 

37 

octagon 

0.18 

0.07 

2.83 

2.6 

Astronomy 

Explorer-17 

188.2 

sphere 

0.621 

0.621 

2 

2 

Scientific 

Sp.  Teles. 

11,000 

cylinder 

112 

14.3 

3.33 

4 

Astronomy 

OSO-7 

634 

9-sided 

1.05 

0.5 

3.67 

2.9 

Solar 

Phys 

OSO-8 

1,063 

cylinder 

5.99 

1.81 

3.76 

4 

Solar 

Phys 

Pegasus-3 

10,500 

cylinder 

264 

14.5 

3.3 

4 

Scientific 

LandSat-1 

891 

cylinder 

10.4 

1.81 

3.4 

4 

Rem  Sens 

ERS-1 

2,160 

box 

45.1 

4 

4 

4 

Rem  Sens 

LDEF-1 

9,695 

1 2-face 

39 

14.3 

2.67 

4 

Environm 

HEAO-2 

3,150 

hexagon 

13.9 

4.52 

2.83 

4 

Astronomy 

Vanguard-2 

9.39 

sphere 

0.2 

0.2 

2 

2 

Scientific 

SkyLab 

76,136 

cylinder 

462 

46.4 

3.5 

4 

Scientific 

Echo-1 

75.3 

sphere 

731 

731 

2 

2 

Comm. 

*Courtesy  of  J.  Wertz  of  Microcosm. 


A.3  Recent  Studies 

According  to  Bowman  and  Moe,19  the  US  Air  Force  has  long  supported  efforts  to  improve  satellite 
orbital  predictions,  in  cooperation  with  the  US  Navy,  NASA,  the  Harvard- Smithsonian  Center  for 
Astrophysics,  and  many  universities  and  Air  Force  contractors.  Over  the  years,  Air  Force  prediction 
efforts  have  evolved  along  two  lines: 


•  In  the  High  Accuracy  Satellite  Drag  Model  (HASDM)20  program,  many  orbiting  satellites  are 
used  to  monitor  the  atmosphere  and  update  an  atmospheric  model  in  real  time 

•  In  the  DMSP  program,  the  spectroscopic  sensors  SSUSI  and  SSULI  will  monitor  the 
atmosphere  and  update  it. 


A  major  objective  of  both  these  programs  is  to  measure  and  predict  absolute  atmospheric  densities. 
Improvement  in  the  knowledge  of  drag  coefficients  contributes  directly  to  that  goal. 

Better  drag  coefficients  also  contribute  to  the  improvement  of  thermospheric  density  models.  A 
landmark  in  the  effort  to  improve  density  models  was  achieved  by  Marcos21  in  1985  when  he 
compared  and  evaluated  14  models  against  accelerometer  measurements  from  seven  satellites.  A 
decade  later,  Chao  et  al.22  used  orbital  data  from  the  ODERACS  I  and  ODERACS  II  families  of 
spheres  to  measure  the  biases  in  the  Jacchia  7 1  and  the  MSIS  90  density  models  near  sunspot 
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minimum.  They  accomplished  this  by  comparing  the  fitted  drag  coefficient,  CD,  with  the  physical 
drag  coefficient,  CDPi  calculated  from  the  actual  momentum  transfer  to  the  satellite.  Results  of  their 
study  indicated  that  the  CDp  values  range  from  2.08  to  2.3  for  altitudes  between  150  km  and  350  km. 
Their  results  also  show  that  MSIS90  model  yields  a  more  accurate  prediction  of  atmospheric  density 
in  that  altitude  region.  More  recently,  Owens23  has  revised  and  improved  the  Jacchia  models  to  create 
the  NASA  Marshall  Engineering  Thermospheric  Model  -Version  2.0  (MET-V  2.0  Model).  This 
model  has  been  subjected  to  an  extensive  statistical  evaluation  by  Wise  et  al.24  Further  recent 
progress  by  Bowman25  in  the  analysis  of  orbital  data  has  refined  the  semiannual  variation  and  other 
parameters  of  the  Jacchia  70  thermospheric  density  model.10  This  new  development  presented  an 
opportunity  to  use  the  improved  Jacchia  model  with  data  from  spherical  satellites  to  refine  our 
knowledge  of  drag  coefficients  and  improve  the  absolute  densities  in  the  model. 

In  2005,  Bowman  and  Moe19  applied  density  corrections  to  the  improved  Jacchia  model  with  tracking 
data  from  numerous  spherical  satellites  to  compute  accurate  fitted  drag  coefficients.  Furthermore, 
they  calculated  physical  drag  coefficients  from  parameters  of  gas-surface  interaction  previously 
measured  in  orbit.  Jacchia  originally  assumed  that  the  physical  drag  coefficient  of  compact  satellites 
is  2.2  at  all  altitudes.  Comparison  of  the  two  types  of  drag  coefficients  has  made  it  possible  to 
remove  from  the  atmospheric  model  the  bias  that  this  assumption  causes. 

Based  on  the  findings  of  Bowman  and  Moe,  the  drag  coefficient,  CD,  of  2.0  should  be  used  with  the 
Jacchia-Bowman  JB2006  model  for  computing  satellite  orbit  decay  and  lifetime.  This  newest  density 
model  considers  the  CD  variation  as  a  function  of  altitude.  For  spacecraft  shape  significantly  different 
from  a  sphere,  such  as  a  cylinder  with  the  relative  wind  along  the  axis  of  rotation,  a  CD  value  of  3 
may  be  used.  A  recent  paper  by  Moe26  gives  the  estimated  ranges  of  CD  for  spheres  and  cylinders  as  a 
function  of  altitude  and  percentage  of  diffuse  reflection. 
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Appendix  B:  25-Year  Lifetime  Predictions  Using  Random  Draw  Approach 


Approaches  have  been  developed  by  the  authors  that  implement  the  solar  and  geomagnetic  modeling 
techniques  presented  in  Section  6.2  of  this  paper.  In  this  first  approach,  the  1  Earth  Research  semi- 
analytic  orbit  propagator  ‘QPROP’  was  used  to  examine  the  7.68  million  cases  contained  in  Table 
B-l.  Each  one  of  these  cases  exhibits  the  type  of  orbital  lifetime  decay  shown  in  Figure  B-l. 


Figure  B-l.  Sample  orbital  decay  profile  of  apogee  and  perigee  altitudes  versus  time. 


For  a  satellite  having  a  ballistic  coefficient  of  200  cm2/kg  and  starting  in  a  circular,  equatorial  orbit  at 
the  altitude  shown,  Figure  B-2  depicts  the  resultant  ranges  of  anticipated  orbit  lifetime.  The 
‘minimum’  and  ‘maximum’  incorporates  the  entire  range  of  orbit  decay  start  times  with  respect  to  the 
solar  cycle  minimum.  The  right-hand  side  of  the  plot  shows  how  variable  the  results  can  get  in  the 
neighborhood  of  25  years’  estimated  lifetime. 
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Orbit  Lifetime  vs  Initial  Circ  Orbit  Altitude 
(CdA/m=200  cmA2/kg;  Equatorial  Orbit) 
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Figure  B-2.  Sample:  orbit  lifetime  (CoA/m  =  200  crrr/kg,  equatorial  orbit)  as  a  function  of  initial  orbit  altitude. 


The  primary  independent  variables  of  the  orbit  lifetime  estimation  process  are  contained  in  Table  B-l. 
The  dependence  of  orbit  lifetime  upon  orbit  inclination  is  shown  for  the  same  200  cm2/kg  sample  case 
in  Figure  B-3.  By  stepping  through  all  of  these  variables  in  the  ranges  and  step  sizes  indicated  in  the 
table,  and  then  detecting  those  cases  that  resulted  in  a  25-year  orbit  lifetime,  the  dependencies 
between  ballistic  coefficient  and  orbit  initial  condition  can  be  found,  as  shown  in  Figures  B-4  and 
B-5.  While  both  the  MSISE2000  and  JB2006  atmosphere  models  are  implemented  in  QPROP,  the 
MSISE2006  model  was  used  for  these  analyses  due  to  its  faster  runtime  with  similar  long-term 
propagation  accuracy.  Random  draws  of  the  triad  of  solar  and  geomagnetic  index  parameters 
(discussed  in  Section  6.2)  were  implemented.  In  order  to  capture  variations  exhibited  by  the  random 
draw  process,  a  number  of  trials  were  used  (four,  in  this  case). 

Table  B-l.  QPROP  Grid  of  Test  Cases 


Parameter 

Lower  Limit 

Upper  Limit 

Step  Size 

Time  into  Solar  Cycle  (days) 

0 

2964.75 

3953/4 

Inclination  (deg) 

0 

90 

30 

CDA/m 

25 

500 

25. 

Perigee  Altitude  (km) 

100 

2000 

50 

Apogee  Altitude  (km) 

250. 

10000 

50 

Number  of  Trials 

0 

3 

1 
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Orbit  Lifetime  Ratio  Between  Equatorial  &  Inclined  Orbits 


Circular  Orbit  Altitude  (km) 

Figure  B-3.  Orbit  lifetime  dependence  upon  orbit  inclination. 

Future  studies  may  use  more  trials  and  incorporate  finer  step  sizes,  but  the  large  computer  runtime 
requirements  of  these  many  cases  led  to  the  initial  selection  of  four  trials  per  initial  set  of  orbit 
conditions.  Through  extensive  simulation,  it  was  found  for  non-sun-synchronous  orbits  that  orbit 
lifetime  results  are  not  sensitive  to  the  three  angular  orbit  elements  (RAAN,  argument  of  perigee,  and 
mean  anomaly)  and  therefore  the  three  initial  values  are  arbitrarily  chosen  and  assumed  for  all  cases. 
Note  that  the  sensitivity  to  RAAN  and  argument  of  perigee  may  be  significant  for  sun-synchronous, 
Magic,  and  Molniya  orbits;  this  is  discussed  in  Section  8.1.  Note:  it  is  recommended  that  sun- 
synchronous  orbit  cases  be  studied  using  a  ‘Method  2’  or  ‘3’  approach  until  such  time  as  their  orbit 
lifetimes  can  be  appropriately  categorized  in  graphical  and/or  functional  form.  Further,  it  was  found 
that  orbits  having  inclinations  greater  than  90  degrees  could  be  well-represented  by  the  pole- 
symmetric  orbits  having  complementary  orbit  inclinations  (justifying  analysis  of  only  0  to  90  degrees 
as  shown  in  Table  B-l). 

The  colored  regions  shown  in  Figures  B-4  and  B-5  denote  the  categorization  of  the  orbit  initial 
conditions  at  the  start  of  the  orbit  decay  with  respect  to  the  IADC  25-year  recommended  post-mission 
lifetime.  The  ‘green’  region  denotes  orbit  initial  conditions  that  will  result  in  an  orbit  lifetime  shorter 
than  25  years  (in  all  observed  cases).  The  ‘yellow’  region  denotes  initial  orbit  conditions  that  could 
result  in  an  orbit  lifetime  that  is  greater  than  the  recommended  25-year  limit  in  certain  circumstances. 

One  can  observe  from  Figures  B-4  and  B-5  that  there  is  a  wide  variety  of  initial  orbit,  timing,  solar 
and  geomagnetic  conditions  that  can  combine  to  produce  an  orbit  lifetime  of  25  years.  These  figures, 
while  helpful,  still  leave  the  user  with  uncertain  knowledge  of  what  the  post-mission  orbit  lifetime 
will  be  specific  to  their  initial  conditions.  Fortunately,  the  results  of  the  7.68  million  analyses  have 
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been  retained;  interpolation  of  these  results  is  possible  to  predict  orbit  lifetime  for  a  specific  set  of 
initial  conditions.  And,  to  the  extent  that  four  sets  of  random  draw  cases  is  not  necessarily  an 
exhaustive  analysis,  additional  cases  can  be  run  to  further  refine  the  orbit  lifetime  estimation  grid  and 
improve  interpolation  results. 


Perigee  vs  Apogee  for  Orbits  Exhibiting  a  25-yr  Lifetime 
(All  orbit  inclinations,  25  <  ballistic  coeff  <  500  cmA2/kg) 
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Figure  B-4.  Perigee  versus  apogee  boundaries  for  25-year  orbit  lifetime  conditions  (25  <  CnA/m  <  500  cm2/kg). 
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Ballistic  Coefficient  vs  Perigee  Altitude  for  25-yr  Orbit  Lifetime 


Figure  B-5.  Ballistic  coefficient  versus  initial  perigee  altitude  for  all  cases  exhibiting  25-year  orbit  lifetime 

(apogee  assumed  <  10,000  km). 

A  by-product  of  the  many  orbit  lifetime  analysis  runs  is  a  set  of  analytical  expressions  (Equation  B-l) 
that  predict  average  orbit  lifetime  as  a  function  of  Hp,  inclination,  ballistic  coefficient,  and  orbit 
eccentricity.  The  term  ‘average’  denotes  orbit  lifetime  averaged  over  all  historical  ranges  of 
solar/geomagnetic  indices  as  well  as  timing  of  the  initial  condition  epoch  with  respect  to  solar 
minimum.  Although  it  is  difficult  to  fit  analytical  expressions  to  an  orbit  lifetime  that  has  been 
averaged  over  all  solar/geomagnetic  historical  conditions  and  over  the  entire  range  of  timings  into  the 
solar  cycle,  a  reasonable  result  has  been  obtained  as  shown  in  Figures  B-6  and  B-7.  Limited  step  size 
resolution  in  Hp,  Ha,  and  time  into  solar  cycle  (four  points  across  cycle)  further  hampered  the  fit 
process. 

This  analytical  approximation  exhibits  a  peak  deviation  of  1 00%  from  the  underlying  estimated 
lifetime  data  in  extreme  cases  (Figure  B-6),  coupled  with  an  average  standard  deviation  of  less  than 
20%  error  above  200  km  (Figure  B-7). 

Orbit  Life(Hp,  Ha,  p,  and  inc)  =  Life(circ  alt  of  Hp)  *  Factor(Havg  ,  inc)  *  Factor(P)  *  Factor(Ecc) 

where: 

Life(circ  alt  of  Hp)  = 

•  if  FTP  <  550  km:  Life(circ  alt  of  Flp)  =  Fi  *  (Ftp  in  km)  2; 

•  else:  Life(circ  alt  of  Flp)  =  Fj  *  (Flp  in  km)  2  *  [  1.0  +  F4  *  (  Flp  -  F3)  5]  (B-l) 
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Factor  (Havg ,  inc)  = 

•  if  Havg  <  155  km:  Factor  (Havg,  inc)  =  1.0; 

•  else  if  Havg  >  550  km:  Factor  (Flavg ,  inc)  =  1.3  -  0.00333333  *  abs  (90  -  inc,  deg); 


•  otherwise:  Factor  (Flavg ,  inc)  = 

-0.9307  *  cos2  (inc)  -  0.2286  *  cos  (inc)  +  2.1011]  exp  [K  *  Havg  ], 


where:  K  =  (2.E-07  *  (90  -  inc,  deg) 2  +  4.E-06  *  abs  (90  -  inc,  deg)  -  0.0009)  if  30 

<  inc  <  150  deg 

or  K  =  -1.33333E-06  *  (inc,  deg)  if  inc  <  30  deg 

or  K  =  -1.33333E-06  *  (180  -  inc,  deg)  if  inc  >  150  deg 


Factor(P)  =  [  F6+ F7  *  (P  in  cm2/kg)  8] 

Factor(Ecc)  =  [  1.0  +  F„  *  (Hp  in  km)2  /  (F10  +  (Hp  in  km)3)  *  (Fu  <F‘2  * ECC>  -  1.0)  ] 


Coeff  Initially,  use: 

Fj  6.350724982689956e-020; 

F2  6.295460674885 14 le+000; 

F3  4.820503248857688e+002; 

F4  3.919897205867653e-006; 

F5  2.3670238903 10601e+000; 

F6  -6.612928883060383e+000; 

F7  2.28569 1478658406e+004; 

F8  -8.49232 1907 172038e-001; 

F9  5.30321 9229 102377e+005; 

Fio  -8.474659230907487e+005; 

F„  1 .00034946260 1627e+000; 

F12  8.79082165932463 le+002; 


If  Resultant  Lifetime  <  1  year: 
6.350724982689956e-020; 
5.946396392270378e+000; 

3. 1853251 3904 1694e+004; 
3.137120106945147e+004; 
3.137356808942190e+004; 

6. 049861 6471 89439e+001; 

1 .203 16106191 7634e+005; 

-9. 1 588280540 1 036  le-001; 
1.522209799457603e+005; 

2.71 807472 1276698e+006; 

1 .00078024 1340044e+000; 
2.065283320597687e+003; 
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Fit  Percent  Error  Versus  Perigee  Altitude 
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Figure  B-6.  Percent  fit  error  versus  perigee  altitude. 


600 


700 


800 


Fit  Statistics  for  Orbit  Lifetime  Fit  Equations 
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Figure  B-7.  Fit  statistics  for  orbit  lifetime  fit  equations. 
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Careful  examination  of  Figures  B-5  and  B-6  shows  that  polar  orbits  experience  reduced  atmospheric 
drag,  due  to  a  polar  orbit’s  reduced  time  spent  under  the  solar  subpoint  (where  solar  heating  causes 
the  atmosphere  to  expand  outward  into  space,  thereby  causing  atmospheric  density  to  rise). 
Therefore,  a  highly-inclined  orbit  requires  a  higher  ballistic  coefficient  (perhaps  an  80  cnr/kg 
increase)  to  match  the  same  25-year  orbit  lifetime  constraint  observed  in  the  low-inclination  orbit 
case  for  an  equivalent  combination  of  perigee  and  apogee  altitudes. 
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Appendix  C:  25-Year  Lifetime  Predictions  Using  Median  Density  Approach 


In  this  approach,  the  satellite  orbit  is  propagated  with  a  semianalytical  tool  (LIFETIME)  using  the 
Jacchia-Bowman  2006  atmosphere  density  model  (JB2006)  based  on  a  50th  percentile  solar  flux 
prediction  (Figure  C-l).  The  new  solar  flux  indices  for  the  JB2006  model  are  also  the  approximated 
predictions  using  the  same  50th  percentile.  The  fundamental  difference  of  this  method  is  that  the  F10.7 
and  Ap  values  are  the  predetermined  13 -month  mean  at  the  50th  percentile,  unlike  the  random  draw 
method  based  on  daily  values  of  several  1 1 -year  cycles.  The  25-year  lifetime  is  determined  by 
iteration  of  initial  eccentricity  holding  initial  apogee  constant.  The  results  of  this  approach  can  be 
used  as  an  independent  verification  of  the  random  draw  approach. 

A  common  epoch  of  1  January  2005  is  assumed  for  all  the  25-year  lifetime  propagations.  This  epoch, 
which  is  near  the  start  of  the  low  solar  flux  period,  gives  the  longest  lifetime.  Thus,  the  results  of  this 
method  are  conservative  and  would  ensure  that  all  cases  did  not  exceed  the  25-year  guideline.  The 
peak-to-valley  variation  in  lifetime  due  to  difference  in  epoch  can  be  as  large  as  4.4  years.  Similar 
sensitivity  to  epoch  is  reported  in  an  early  study  by  Meyer  and  Chao.27 


Year  from  1  Jan  2007 


Figure  C-l.  A  50th  percentile  of  13-month  mean  solar  flux  (F10.7). 

A  large  number  of  25-year  cases  were  generated  based  on  the  following  assumed  orbit  conditions. 
Again,  it  was  found  for  non-sun-synchronous  orbits  that  orbit  lifetime  results  are  not  sensitive  to  the 
three  angular  orbit  elements  and  therefore  the  three  initial  values  are  arbitrarily  chosen  and  assumed 
for  all  cases.  Note  that  the  sensitivity  to  RAAN  and  argument  of  perigee  may  be  significant  for  sun- 
synchronous,  Magic  and  Molniya  orbits;  this  is  discussed  in  Section  8.1. 

•  Apogee  altitude  =  500  km  to  1200  km  at  100  km  step 

•  Ballistic  coefficient  =  100,  300,  500  cm2/kg 

•  Inclination  =  0,  30,  60,  90  deg 

•  Right  ascension  of  ascending  node  =  260  deg 

•  Argument  of  perigee  =  80  deg 

•  Mean  anomaly  =  0  deg 
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Figures  C-2  and  C-3  show  the  converged  solutions  for  25-year  lifetimes  at  90  and  30  deg  inclinations, 
respectively.  The  tolerance  for  a  converged  solution  is  when  the  lifetime  is  within  +/-  6  months  of  25 
years.  It  is  interesting  to  see  that  the  three  curves  for  the  three  values  of  ballistic  coefficient  are  nearly 
linear  for  apogee  altitude  higher  than  600  km.  The  well-behaved  curves  in  Figures  C-2  and  C-3 
suggest  that  the  solutions  may  be  more  conveniently  represented  by  families  of  semiempirical 
equations.  A  general  equation  with  empirical  constants  has  been  derived  as: 

Eccentricity  =  Ci(Ha  -  Bi)  -  C2/(Ha  -  B2)  (C-l) 

where  Ha  is  initial  apogee  altitude  and  the  four  constants,  Ci  Bi  C2  and  B2,  are  determined  by  fitting  to 
each  curve  in  the  two  figures.  The  first  term  in  the  above  equation  gives  the  linear  change,  while  the 
second  term  approximates  the  gradual  deviations  at  low  apogee  region.  The  values  of  these  constants 
for  all  the  above  cases  are  given  in  Table  C-l. 


Initial  apogee  altitude,  km 


Figure  C-2.  Eccentricity  vs.  initial  apogee  altitude  for  25-year  lifetime  at  90  deg  inclination 
with  CDA/m  =  100,  300  and  500  cnr/kg. 
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Figure  C-3.  Eccentricity  vs.  initial  apogee  altitude  for  25-life  lifetime  at  30  deg  inclination 
with  CnA/m  =  100,  300  and  500  cm2/kg. 

Table  C-l.  Values  of  Ci,  Bi,  C2  and  B2  for  Equation  4 


Inclination, 

deg 

CD  A/m,  cnf/kg 

Ci 

Bi 

c2 

b2 

0 

100 

0.00008 

480 

0.4 

540 

300 

same 

550 

0.65 

600 

500 

same 

580 

0.6 

640 

30 

100 

0.00008 

465 

0.42 

550 

300 

same 

550 

0.46 

608 

500 

same 

585 

0.45 

650 

60 

100 

0.00008 

450 

0.75 

470 

300 

same 

515 

0.95 

545 

500 

same 

560 

0.87 

595 

90 

100 

0.00008 

435 

0.76 

500 

300 

same 

520 

0.87 

540 

500 

same 

550 

0.8 

600 

This  method  with  the  semi-empirical  equation  (Eq.  C-l)  provides  a  convenient  way  to  determine  the 
initial  eccentricity  for  achieving  a  25-year  lifetime  based  on  orbit  and  spacecraft  conditions  at  end-of- 
life.  The  required  AV  or  propellant  for  orbit  insertion  can  then  be  computed  (see  Figure  C-4  as  an 
example).  The  value  of  eccentricity  determined  by  this  method  may  be  verified  using  either  Method  1 
(previously-generated  results)  or  Method  2  (semi-analytical  integration)  discussed  in  this  report. 
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Figure  C-4.  AV  vs.  initial  apogee  altitude  for  achieving  25-year  lifetime  at  30  deg  inclination 
with  CDA/m  =  100,  300  and  500  cm2/kg. 
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Appendix  D:  Orbit  Lifetime  Sensitivity  to  Sun-Synchronous,  Magic,  Molniya 
Orbits  and  Thermospheric  Global  Cooling 


D.1  Orbit  Lifetime 

For  sun-synchronous  orbits,  orbit  lifetime  has  some  sensitivity  to  the  initial  value  of  RAAN  due  to  the 
density  variations  with  the  local  sun  angle.  Results  of  25-year  propagations  using  both  JB2006  and 
MSIS2000  show  similar  trends — that  orbits  with  6:00  a.m.  local  time  have  longer  lifetime  than  orbits 
with  12:00  noon  local  time  by  about  500  days.  This  maximum  difference  (500  days)  translates  into  a 
5%  error.  It  is  recommended  that  the  long-term  (25  year)  lifetime  prediction  of  a  sun-synchronous 
orbit  should  be  performed  only  by  a  semi-analytical  or  numerical  integration  tool. 

For  Molniya  orbits,  iteration  has  been  found  to  be  difficult  when  computing  the  25-year  lifetime  due 
to  the  coupling  in  eccentricity  between  the  third-body  perturbations  and  the  drag  decay.  Figure  D-l 
shows  the  15-year  history  of  perigee  altitude  of  a  particular  Molniya  orbit  with  an  initial  RAAN  of  97 
deg.  The  long-term  (7-year)  sinusoidal  variation  is  caused  by  the  third-body  attractions  on  orbit 
eccentricity.  The  perigee  variation  can  be  out  of  phase  by  90  degrees  if  the  initial  RAAN  becomes  7 
deg.  The  atmospheric  drag  force  becomes  significant  only  when  the  perigee  is  considerably  below 
800  km.  The  strong  sensitivity  to  initial  RAAN  and  epoch  makes  it  difficult  to  estimate  the  25 -year 
lifetime  using  empirical  equations.  Similar  to  the  Sun-synchronous  orbits,  it  is  recommended  that 
Method  2  or  3  be  used  when  assessing  the  initial  conditions  that  will  achieve  a  25 -year  lifetime  for 
this  type  of  orbit. 


Years  after  epoch 


Figure  D-l.  Perigee  altitude  history  of  a  Molniya  orbit. 

The  so-called  ‘Magic’  orbit  is  defined  as  having  a  three-hour  period  with  a  repeating  ground  track, 
where  the  orbit  is  highly-elliptical,  critically-inclined  and  Sun-synchronous.  A  typical  Magic  orbit  is 
given  below: 

a  =  10572.87  km 
e  =  0.35 
i  =  1 16.565  deg 
D.  =  45  deg 
ca  =  270  deg 
M  =  0  deg 
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Based  on  the  above  orbit,  the  maximum  perigee  altitudes  (for  not  exceeding  25-year  lifetime)  for 
three  selected  ballistic  coefficients  (100,  300,  500  cm2/kg)  were  determined  using  LIFETIME  and 
plotted  against  local  mean  solar  time  of  ascending  node  in  Figure  D-2.  The  strong  sensitivity  to  initial 
RAAN  is  due  to  the  sun-synchronous  nature  of  the  Magic  orbit.  The  difference  in  initial  perigee 
altitude  of  the  disposal  orbits  can  be  as  large  as  180  km,  depending  on  the  difference  in  local  mean 
solar  time. 


550 
500 
450 
400 
350 
300 
250 
200 

0  2  4  6  8  10  12  14  16  18  20  22  24 
Local  mean  solar  time,  hours 
Figure  D-2.  Estimated  perigee  altitudes  vs.  local  mean  solar  time  of  ascending  node  for  Magic  orbits. 

Although  the  above  figure  (Figure  D-2)  can  serve  as  a  guide  to  determine  the  perigee  altitude  of  a 
disposal  orbit  of  a  Magic  orbit,  it  is  again  recommended  that  Method  2  or  3  be  used  when  determining 
the  25 -year  lifetime  for  this  type  of  orbit. 

D.2  Implication  of  Observed  Thermospheric  Global  Cooling 

Recent  indications  of  global  cooling  in  the  thermosphere  may  have  a  gradually  increasing  role  in  orbit 
lifetime  estimation.  The  thermosphere  is  defined  to  occur  roughly  between  80  and  500  km  altitude, 
which  is  a  key  part  of  the  LEO  regime  for  which  the  ISO  standard  is  being  developed.  Both  satellite 
measurements'  and  theoretical  models  ’  indicate  that  the  thermosphere  is  cooling  off,  causing 
density  to  lower.  The  mechanism  causing  this  change  appears  to  be  that  as  C02  concentrations  have 
increased  (from  320  ppmv  in  1965  to  around  380  ppmv  in  2005)  at  altitudes  below  30km,29  the  upper 
atmosphere  is  correspondingly  cooling  down.  It  is  estimated  that  because  of  this  effect,  atmospheric 
density  will  decrease  by  between  1.7%  and  2%  ’  per  decade. 

To  study  the  implications  of  this  finding,  we  placed  a  1.7%  per  decade  decrease  into  the  LIFETIME 
and  QPROP  software  programs.  The  results  are  shown  in  Figures  D-3  and  D-4.  In  Figure  D-3,  the 
orbit  decay  profiles  for  the  equatorial  (0  degree  inclination)  and  polar  (90  degree  inclination)  cases 
are  compared  for  the  current  atmosphere  model  predictions  (profiles  ‘A’)  and  the  atmosphere  models 
with  an  overlaid  thermospheric  cooling  density  effect  of  1.7%  per  decade.  Figure  D-4  contrasts  with 
the  sample  case  used  in  Figure  B-2.  In  both  cases,  note  that  the  introduction  of  the  1.7%  per  decade 
thermospheric  cooling  effect  density  decrease  led  to  an  orbit  lifetime  increase.  Specific  to  the  IADC- 
recommended  25-year  orbit  lifetime  criterion,  Figure  D-3  shows  a  predicted  lifetime  increase  of 
approximately  6-7%. 
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Ratio  of  Orbit  Lifetimes  Apogee  Altitude,  km 
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Figure  D-3.  Orbit  lifetime  estimates  including  and  excluding  thermospheric  cooling 
for  0°  &  90°  inclination  cases. 
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Figure  D-4.  Orbit  lifetime  variation  induced  by  thermospheric  global  cooling  effect 
on  sample  case  from  Fig.  B-2. 


Orbit  Lifetime  Ratio  of  Thermospheric  Cooling  vs  non-TC 
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Appendix  E:  Discussion  of  Conditional  Probability 


The  reliability  requirement  given  in  6.2,  “The  space  system  shall  be  designed  such  that  the  probability 
of  successful  end-of-mission  disposal,  including  depletion  of  energy  sources,  equals  or  exceeds  0,9  at 
the  time  disposal  is  executed,”  represents  a  typical  and  consistent  approach  to  managing  the 
multiphased  mission  from  a  probability  point  of  view. 

Assuming  the  orbital  scenario  is  constituted  by  the  following  two  main  phases: 

•  Nominal  mission  (named  M),  and 

•  Disposal  (named  D). 


the  Conditional  Probability  Theorem  allows  us  to  state 

P(M  nD)  =  P(M)xP(D/M)  (E-l) 

It  follows  that 

p{D/MhPimf  ,E'2) 

Therefore,  the  probability,  P{D/M),  of  successfully  performing  the  end-of-mission  disposal  once  the 

nominal  mission  has  been  completed  may  be  estimated  starting  from: 

•  P(M)  =  mission  success  probability  (normally  evaluated  in  the  frame  of  the  dependability 
program  to  verify  the  compliance  of  the  developed  design  with  regard  to  the  applicable 
mission  reliability  requirement),  and 

•  P(MnD)  =  probability  of  correctly  performing  both  mission  and  disposal  phases.  (This  is 
generally  an  extension  of  the  P{M)  model.) 

The  following  cases  are  presented  to  clarify  the  application  of  Equation  (E-2). 

E.1  Case  1 


The  end-of-mission  disposal  is  performed  making  use  of  “all”  subsystems  involved  in  the  nominal 
mission,  as  illustrated  in  the  Venn  diagram  shown  in  Figure  E-l. 
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Figure  E-l.  All  subsystems  involved  in  the  nominal  mission  required  for  disposal. 


In  this  way,  the  reliability  model  defined  for  the  nominal  mission  [0;  Tmissio„\  is  extended  to  cover  both 
phases  [0,  T mission  T disposal}  • 


P(D/M)  = 


R 


system 


(T, 


+  T 

mission  disposal  > 


R 


(T.) 

\  mission  / 


system  V  mission 


(E-3) 


E.2  Case  2 

The  end-of-mission  disposal  is  performed  making  use  of  “some”  subsystems  involved  in  the  nominal 
mission,  as  illustrated  in  the  Venn  diagram  shown  in  Figure  E-2. 


Figure  E-2.  Some  subsystems  involved  in  the  nominal  mission  are  used  for  disposal. 


In  this  case,  the  function  R  ’ svs,em(j)  includes  the  reliability  of  the  subsystems  used  for  disposal 
estimated  over  the  interval  [0;  Tmission  +  Tdisposat\.  For  the  remaining  subsystems,  the  interval  [0;  Tmission] 
is  considered.  For  all  subsystems,  the  reliability  models  defined  for  mission  success  may  be  used, 
even  if  with  different  time  durations. 


The  reliability  function  i?system(t)  refers  to  the  nominal  mission  (as  presented  in  Case  1). 
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E.3  Case  3 


The  end-of-mission  disposal  is  performed  making  use  of  dedicated  equipment  and  “some” 
subsystems  involved  in  the  nominal  mission,  as  illustrated  in  the  Venn  diagram  shown  in  Figure  E-3. 


Figure  E-3.  Disposal  requires  some  systems  used  for  normal  operations  plus  some  systems 

specifically  used  for  disposal  only. 


In  this  extension  of  Case  2,  the  R"system(t )  function  includes  the  contribution  of  the  equipment  for 
disposal:  considered  in  active  or  passive  mode  (as  applicable)  in  the  interval  [0;  Tmissio„\  and  in  active 
mode  in  the  interval  [Tmi-55I-on,  TmiSSion  T  T dispose*  t\  • 

The  reliability  function,  Rsystem(t),  refers  to  the  nominal  mission  (as  presented  in  Case  1). 
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Several  cases  may  be  described.  The  cases  presented  in  this  appendix  are  provided  to  clarify  the 
approach  taken  in  this  standard. 
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2.  Send  to  the  Preparing  Activity  specified  in  block  8. 

NOTE:  Do  not  use  this  form  to  request  copies  of  documents,  or  to  request  waivers,  or  clarification  of  requirements  on 
current  contracts.  Comments  submitted  on  this  form  do  not  constitute  or  imply  authorization  to  waive  any  portion  of 
the  referenced  document(s)  or  to  amend  contractual  requirements.  Comments  submitted  on  this  form  do  not 
constitute  a  commitment  by  the  Preparing  Activity  to  implement  the  suggestion;  the  Preparing  Authority  will 
coordinate  a  review  of  the  comment  and  provide  disposition  to  the  comment  submitter  specified  in  Block  6. 
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